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Abstract 

We present measurements of decay amplitudes and triple-product correlations in B —* <f>K* 
decays based on 253 fb _1 of data recorded at the T(45) resonance with the Belle detector at 
the KEKB e + e _ storage ring. The decay amplitudes for the three different helicity states are 
determined from the angular distributions of final state particles. The longitudinal polarization 
amplitudes are found to be 0.45 ±0.05 ±0.02 for B° <f)K*° and 0.52 ±0.08±0.03 for B + —> <j)K* + 
decays. CP- and T-odd CP-violating triple-product asymmetries are measured to be consistent 
with zero. 

PACS numbers: 13.25.Hw, 11.30.Er 
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The vector-vector B —> <pK* decay processes provide clear insights into the underlying 
b —> s transition by virtue of their clear experimental signatures and relatively unambiguous 
theoretical interpretation, especially of the many angular correlations that can be formed 
among the final-state particles. The decays are described by second order penguin diagrams, 
the first order b —> s transition being forbidden in the Standard Model (SM). The angular 
information allows the CP-even and CP-odd states that comprise the B° —> decay 

to be distinguished. Our previous measurement jjj and a recent report by BaBar Q both 
suggest that the longitudinal polarization component differs from predictions based on the 
factorization assumption. 

In this letter we report on a further study of this anomaly that is based on a larger data 
sample and uses observables that are expected to be sensitive to the effects of new physics. 
We present the first full three-dimensional angular analysis for B + —> <pK* + and an extended 
study for B° —■> (f>K*°. The decay modes 0 —■> K + K ~, K*° —■> K + ir~, K* + —■> K° s tt + , 
and K* + —> K +/ ii 0 are considered. Charge conjugate modes are implied everywhere unless 
otherwise specified. We report measurements of direct CP asymmetries, triple-product 
correlations and related T-odd CP -violating asymmetries !3], and other observables that 
are sensitive to New Physics (NP) 4i], 

This analysis uses a data sample that contains 275 x 10 6 BB pairs collected on the T(4S') 
resonance by the Belle detector I at the KEKB e + e collider jf|. The Belle detector is 
a general purpose magnetic spectrometer equipped with a 1.5 T superconducting solenoid 
magnet. Charged tracks are reconstructed in a central drift chamber (CDC) and a silicon 
vertex detector (SVD). Photons and electrons are identified using a CsI(Tl) electromagnetic 
calorimeter (ECL) located inside the magnet coil. Charged particles are identified using 
specific ionization ( dE/dx ) measurements in the CDC as well as information from aerogel 
Cherenkov counters (ACC) and time of flight counters (TOF). 

Event reconstruction is performed as described in Ref. |l|. Candidate B mesons are 
reconstructed from (j) and K* candidates and are identified by the energy difference A E = 
£™ s — FbSfm, beam constrained mass M hc = yj (FbSm) 2 ~~ ( P f 18 ) 2 , and K + K~ invariant 
mass ( M k + k -), where E ^ m is the beam energy in the center-of-mass system (cms), and 
F|j ms and are the cms energy and momentum of the reconstructed B candidate. The B- 
rneson signal region is defined as M bc > 5.27 GeV/c 2 , |AP| < 45 MeV, and \M K + K - — M^\ < 
10 MeV/c 2 . The invariant mass of the K* —> Kit candidate is required to be less than 70 
MeV/c 2 from the nominal K* mass. The signal region is enlarged to —100 MeV < A E < 
80 MeV for B + —> <f)K* + (K* + —> K + tt °) because of the effects of shower leakage on the A E 
resolution. An additional requirement cos 6k* < 0.8 is applied to reduce low momentum 
7T° background, where 6k* is the angle between the direction opposite to the B and the 
daughter kaon in the rest frame of K*. These requirements do not effect our results based 
on a MC study. In the signal region, about 1% of the events have multiple candidates. The 
candidate with the smallest y 2 value from B vertex finding and the best 7T° mass in the 
Ii* + —> A' + 7t° decay is used. 

The dominant background is e + e _ — > qq (q = u, cl, c, s ) continuum production. Several 
variables including S± [7|, the thrust angle, and the modified Fox-Wolfram moments defined 
in Ref. Q are used to exploit the differences between the event shapes for continuum qq 
production (jet-like) and for B decay (spherical) in the cms frame of the T(4S'). These 
variables are combined into a single likelihood ratio 1Z S = C S /(C S + C q q), where C s ( C qq ) 
denotes the signal (continuum) likelihood. The selection requirements on 1Z S are determined 
by maximizing the value of N S /^N S + Nj, in each R-flavor-tagging quality region [9j|, where 
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N s (Nb) represents the expected number of signal (background) events in the signal region. 

Backgrounds from other B decay modes such as B —> K + K~K*, B —> /o(980)A'*(/ 0 —>► 
Ii + K~), B —> cpK 7r, A —> K + K~Ktt, and cross-feed between the 0A'* and <pK decay 
channels are studied. The contributions from B —> K + K~K* and B —> foK* are estimated 
from a fit to AA, Mb c , and Mk+k - distributions. The M k +k- distribution for B —> 
K + K~K* is determined from Monte Carlo (MC) simulations. The / 0 (980) line shape is 
obtained from MC, where an A-wave Breit-Wigner with a 61 MeV/c 2 intrinsic width jlCj] is 
assumed. The uncertainty in the / 0 (980) width (40-100 MeV/c 2 ) H3 is taken as a source 
of systematic error. The contributions from B —* K + K~K* ( B —> foK*) are estimated 
together with the <pK* signal and are found to be 1 to 7% (1 to 3%) [l2j of the signal yield. 
The background from B —> (pKn decays is evaluated with fits to the Kn invariant mass 
and is found to be about 1%. The contamination from four-body B —> K + K~Kn decays is 
checked by performing fits to the events in the <p —> K + K~ and K* —> Kti mass sidebands 
and is found to be less than 1%. To remove the contamination from <pK decays, these decays 
are explicitly reconstructed and rejected. 

The signal yields (N s ) are extracted by extended unbinned maximum-likelihood fits per¬ 
formed simultaneously to the A E, M^ c and Mk+k- distributions. Reconstructed B candi¬ 
dates with |AA| < 0.25 GeV, M^ c > 5.2 GeV/c 2 , and Mk+k- < 1-07 GeV/c 2 are included 
in the fits. The signal probability density functions (PDFs) are a single Gaussian in Mb c , 
a core Gaussian plus a Bifurcated Gaussian (Gaussian with different widths on either side 
of the mean) as the tail in the AA distribution, and a Breit-Wigner shape in Mk+k-- The 
means and widths of A A and M bc are verified using B —> J/ ip K* decays. The mean and 
width of the <p mass peak are determined using an inclusive (p -+ K + K~ data sample. 

The PDF shapes for the continuum events are parameterized by an ARGUS function 
[13] in Mbc, a linear function in AA, and a sum of a threshold function and a Breit-Wigner 
function in Mk+k~ ■ The parameters of the functions are determined by a fit to the events 
in the sideband. The signal and background yields are allowed to float in the fit while other 
PDF parameters are fixed. The direct CP asymmetries, A C p = > are also 

l+yB—rj J+iV (B—*j ) 

studied. The measured signal yields and direct CP asymmetries are summarized in Table HI 
The distributions of AA, Mb c , and M k +k- are shown in Fig. |T] 
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FIG. 1: Distributions of the AA, Mb c , and M k +k~ f° r B° (pK*° ((a),(b) and (c)), and 
for B + -+ cpK* + ((d),(e) and (f)), with other variables in signal region. Solid curves show the fit 
results. The continuum background components are shown by the dashed curves. The dark shaded 
areas represent the contributions from B —+ K + K~K* and B -+ foK* decays. 
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TABLE I: Number of events observed in the signal region (iV ev ), signal yields (N s ) and the direct 
CP asymmetries ( Acp ) obtained in the fits, with statistical and systematic uncertainties. 


Mode 

N ev N s Acp 

0A'*° 

309 173 ± 16 0.02 ± 0.09 ± 0.02 

0A'*+ 

K* + (Kg7T + ) 

K* + (K+tt°) 

173 85l^ -0.02 ±0.14 ±0.03 

76 37.91^0 -0.14 ±0.21 ±0.04 
97 47.3l|i 0.09 ±0.19 ±0.04 


The decay angles of a 5-meson decaying to two vector mesons 0 and K* are defined in 
the transversity basis 3|- The x-y plane is defined to be the decay plane of K* and the x 


axis is in the direction of the 0-meson. The y axis is perpendicular to the x axis in the decay 
plane and is on the same side as the kaon from the K* decay. The z axis is perpendicular to 
the x-y plane according to the right-hand rule, 9 tI (0tr) is the polar (azimuthal) angle with 
respect to the z-axis of the K + from 0 decay in the 0 rest frame, and 9k* is defined earlier. 

The distribution of the angles, 9k*, 9 tr , and 0 tr is given by Til 


(0tr, COS 6> tr , cos 9 k * ) _ 9 
C?0trd COS 9t r d COS 9k* 327T 

|A_l| 2 2 cos 2 9 tI sin 2 9 K * 

+ |A|| | 2 2 sin 2 9 tI sin 2 0 tr sin 2 9 K * 

+ | A)| 2 4 sin 2 9 tI cos 2 0 tr cos 2 9 K * 

± v^Re^Ao) sin 2 #tr sin 20 tr sin 2 9k* 

— yV^lm^A^Ap) sin 29 tr cos 0 tr sin 29k* 

— 2?/Im(A|Aj_) sin 2 9 tr sin 0 tr sin 2 9 K *} , (1) 

where A$, An, and A± are the complex amplitudes of the three hclicity states in the transver¬ 
sity basis with the normalization condition | Ao | 2 + | A ||| 2 + |71_l| 2 = 1 , and y = +1 (— 1 ) 
corresponds to B ( B ) mesons and is determined from the charge of the kaon or pion in the 
K* decay. The longitudinal polarization component is denoted by A 0 ; A± (Ay) is the trans¬ 
verse polarization along the z-axis (y-axis). The value of |A_|_| 2 (|A 0 | 2 + l^||| 2 ) is th e CP-odd 
(CP-even) fraction in the decay B —* 0/i*° ^|. The presence of final state interactions 
(FSI) results in phases that differ from either 0 or ± n. 

The complex amplitudes are determined by performing an unbinned maximum likelihood 
fit to the B —> (j>K* candidates in the signal region. The combined likelihood is given by 


£ = n e (^tr, cos 9 tr , cos 9 k *) Y, fj R j (0tr, cos 0 tr , cos 9 K * ) , 
i j 


( 2 ) 


where j denotes the contributions from 0/i*, qq, K + K~K* and foK*', Rj is the angular 
distribution function (ADF). The ADF R q g is determined from sideband data, and Rk+k-k* 
from events with 1.04 GeV/c 2 < M k +k- < 1-075 GeV/c 2 ; Rj 0 k* is obtained from B —> f$K* 
MC events. The detection efficiency (e) is determined using MC simulations assuming 
a phase space decay. The fractions fj are parameterized as a function of A E, M hc and 
M k +k~■ The value of arg(A 0 ) is set to zero and |A||| 2 is calculated from the normalization 
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condition. The four parameters (|A 0 | 2 , |Aj_| 2 , arg(Ay), and arg(AjJ) are determined from 
the fit. There is a two-fold ambiguity in the solutions for the phases; the chosen set of 
solutions is the one suggested in Ref. 0 Figure El shows the angular distributions with 
projections of the fit superimposed. The obtained amplitudes are summarized in Table ITU 
The systematic uncertainties on the amplitudes are dominated by the efficiency modeling 
(4-5%), continuum background (3-4%), slow pion efficiency (2-3%), and K + K~K* ADF 
(1-2%). The remaining possible systematic errors, such as the angular resolution, signal 
yields, background from higher K* states, and width of the / 0 , are estimated to be less than 
1 %. 


TABLE II: The decay amplitudes obtained for B° —> cj)K*° and B + —► cpK* + . The first uncertain¬ 
ties are statistical and the second are systematic. 


Mode 

O 

* 

4>K*+ 

|A 0 | 2 

0.45 ± 0.05 ± 0.02 

0.52 ±0.08 ±0.03 

A±| 2 

0.30 ± 0.06 ± 0.02 

0.19 ±0.08 ±0.02 

arg(An) (rad) 

2.39 ±0.24 ±0.04 

2.10 ±0.28 ±0.04 

arg(Aj) (rad) 

2.51 ±0.23 ±0.04 

2.31 ±0.30 ±0.07 



FIG. 2: Projected distributions of the three transversity angles for B° —> <j)K*° ((a),(b) and (c)), 
and for B + —> (j)K* + ((d),(e) and (f)). Solid lines show the fit results. The points with error bars 
show the efficiency corrected data after background subtraction. The two I\* + decay modes are 
combined in (d),(e) and (f). The discontinuity in (d) is due to the requirement of cos 6%* < 0.8 in 
B+ -»• (j)K* + (K* + -»• K+tt°). 

The triple-product for a B meson decay to two vector mesons takes the form if • (ei x ej), 
where if is the momentum of one of the vector mesons, and e? and ef are the polarizations 
of the two vector mesons. The following two T-odd |3| |l7[ quantities 

= Im(A_i_Ag) , = Im(A_i_Ap , (3) 

provide information on the asymmetry of the triple products. The SM predicts very small 
values for Ay and Ay. The comparison of these triple product asymmetries (Ay and Ay) 
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with the corresponding quantities for the CP-conjugate decays (A^ and A provides an 
observable sensitive to T-odd CP-violation. 

Additional variables that can be accessed by angular analyses are suggested in Ref. 4] 
and are given by 


A_u = -Im (A ± A* - Aj_A*), A No = Re(A||A* + A\\A* 0 ), 

= -lm(A ± A* + A ± A*), E,| 0 = Re^A*-A l{ A* 0 ), 

Aaa = |(|A A | 2 + p A | 2 ), E AA = i(|A A | 2 - q A | 2 ), (4) 

where the subscript A is either 0, ||, or A, and i is 0 or ||. The variables Aj_ 0 and Aj_j,| are 
sensitive to T-odd CP -violating new physics. The following equations should hold in the 
absence of NP: 


Saa = 0 , E|| 0 = 0 , A ±J = 0 . (5) 

By separating B° and B° samples and rearranging fitting parameters in the unbinned 
maximum likelihood fit, we measured the decay amplitudes for the B° and P°, the triple¬ 
product correlations, and the other NP-sensitive observables as given in Table II I II and IIVI 
The T-odd CP-violating variables A^ 0 and Aj_|| are measured to be 0.± 0.03 and 
0.01 ± 0.10 ± 0.02, respectively, consistent with the SM predictions. 


TABLE III: The measured decay amplitudes and triple-product correlations in the B° and B° 
samples. 


Mode 

B° 



B° 



IA)I 2 


0.39 ±0.08 ±0.03 

0.51 

± 

0.07 

± 

0.02 

\A±\ 2 


0.37 ± 0.09 ±0.02 

0.25 

± 

0.07 

± 

0.01 

arg(A||) (rad) 

2.72^38 ± 0.14 

2.08 

± 

0.31 

± 

0.04 

arg(A^) (rad) 

2.81 ±0.36 ±0.11 

2.22 

± 

0.35 

± 

0.05 

A° T 


0.13loa4 ± 0.04 

0.28 

± 

0.08 

± 

0.01 



0.03 ±0.08 ±0.01 

0.03 

± 

0.06 

± 

0.01 


In summary, improved measurements of the decay amplitudes for B —> <f>K *, based on 
fits to angular distributions in the transversity basis, are presented. The results are con¬ 
sistent with our previous measurements jH with improved precision. The measured value 
of |A_l | 2 shows that CP-odd (|Aj_| J ) and CP-even (|A 0 | 2 + |A||| 2 ) components are present 
in (j)K * decays in a ratio of about 1:2. Phases of both Aj_ and A\\ differ from zero or —i r 
by 4.3 standard deviations (cr), which provides evidence for the presence of final state in¬ 
teractions. The measured direct CP asymmetries in these modes are consistent with zero; 
the corresponding 90% confidence level limits are —0.14 < Acp(4>K* 0 (K + 7i ~)) < 0.17, and 
—0.25 < Acp(4>K * + ) < 0.22. Measurements of the T-odd CP-violation sensitive differences 

between triple product asymmetries, A ^ and A^ — indicate no significant devia¬ 

tions from zero, consistent with BaBar measurements fl]. Our data shows no significant 
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TABLE IV: A and £ values obtained from the decay amplitudes measured for B° and B° separately. 


Aoo = 

0.45 

± 

0.05 

± 

0.02 

Soo — 

-0.06 

± 

0.05 

± 

0.01 

A|IH = 

0.24 

± 

0.06 

± 

0.02 

£|IH = 

-0.01 

± 

0.06 

± 

0.01 

A_l_l = 

0.31 

± 

0.06 

± 

0.01 

£j__l = 

0.06 

± 

0.05 

± 

0.01 

£±o = 

- 0 . 

41 

■fO.16 

-0.14 

± 

0.04 

Aj_o = 

0 . 

16 

1-0.16 

-0.14 

± 

0.03 

S UI = 

-0.06 

± 

0.10 

± 

0.01 

A ±|| = 

0.01 

± 

0.10 

± 

0.02 

A l|o = 

-0.45 

± 

0.11 

± 

0.01 

S l|o 

-0.11 

± 

0.11 

± 

0.02 


deviations from the expectations: £aa = 0, Eyo = 0, and Aj_j = 0, indicating no evidence 
for new physics. 
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